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The electric dipole moment (EDM) of a particle violates both time reversal (T) and space 
reflection (P) symmetries. There have been recent suggestions for searches of the electron EDM 
using solid state experiments [[jj,^ . These experiments could improve the sensitivity compared 
to present atomic and molecular experiments by several orders of magnitude. In the present 
paper we calculate the expected effect. We also suggest that this kind of experiment is sensitive 
to T,P- violation in nuclear forces and calculate effects caused by the nuclear Schiff moment. 

The compounds under consideration contain magnetic Gd 3+ ions and oxygen O 2- ions. We 
demonstrate that the main mechanism for the T,P-odd effects is related to the penetration of 
the Oxygen 2p-electrons to the Gd core. All the effects are related to the deformation of the 
crystal lattice. 
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I. INTRODUCTION 

Violation of the combined symmetry of charge conjugation (C) and parity (P) has been discovered in the decay 
of the K° meson about 40 years ago || . The exact origin of this symmetry violation remains an enigma, although 
the so called standard model of electroweak interactions can describe these processes phenomenologically. It 
has also been proposed by Sakharov Q that the matter-antimatter asymmetry observed in our Universe could 
have arisen from a CP-violating interaction active at an early stage of the Big Bang. The CP-violation implies 
a time-reversal (T) asymmetry and hence violation of both T- and P-symmetries, because there are strong 
reasons to believe that the combined CPT-symmetry should not be violated || . Electric dipole moment (EDM) 
of a system in a stationary quantum state indicates violation of T- and P-symmetries. This is why searches for 
EDM of elementary particles, atoms and molecules are very important for studies of violations of fundamental 
symmetries |J. A property that is closely related to EDM is the so called Schiff moment, see e.g. Ref. ||. A 
non zero Schiff moment also indicates T- and P-violation. In the present paper we consider T,P-odd effects in 
solids due to the EDM of the electron and due to the nuclear Schiff moment. 

The present best limitation on the electron EDM comes from the Berkeley experiment with an atomic Thallium 
beam 0, 

d e < 1.6 x l(T 27 e cm. (1) 

There are some ideas for improvement of the sensitivity. One way of improvement is an experiment with 
metastable levels of PbO molecules |J. A breakthrough could be achieved in solid state experiments with 
compounds containing uncompensated spins. This idea was suggested by Shapiro in 1968 ||. Application of 
a strong electric field to electrons bound within a solid would align the EDMs of these electrons. This should 
lead to a simultaneous alignment of the electron spins; the magnetic field arising from this alignment could 
be detected experimentally. Another possibility is to polarize electrons by the external magnetic field. This 
causes alignment of electron EDMs, and hence induces a voltage across the sample that could be detected. 
An experiment of this kind has been performed with nickel-zinc ferrite |ni(| . However due to experimental 
limitations the result was not very impressive. Interest in this approach has been renewed recently due to a 
suggestion of Lamoreaux |IJ and Hunter j|] to perform similar experiments with Gadolinium Gallium Garnet, 
GdaGa50i2 (GGG), and Gadolinium Iron Garnet GdsFesO]^ (GdIG) employing new experimental techniques. 
The estimates of sensitivity look highly promising. In the present work we calculate the expected effects. 
The best limitation on the Nuclear Schiff moment (NSM) comes from the Seattle experiment with atomic 
Hg @, 
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S N ( 199 Hg) < 0.7 x 10- 49 e cm 3 . (2) 

The expected value of Sn expressed in terms of the fundamental CP-violating interaction depends on the 
structure of the particular nucleus, however for heavy nuclei this dependence is not that strong and therefore 
we will take the value (|^) as a reference point. Later, we will comment on a possible additional enhancement 
due to a special structure of the 155 Gd nucleus. One can search for NSM in the experiments with GGG and 
GdIG simultaneously with searches for the electron EDM. The point is that due to the hyperfine interaction 
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nuclear polarization causes polarization of electrons and vice versa. Another possibility is to use a compound 
without unpaired electrons, for example Lutetium Gallium Garnet, Lu3GasOi2 or Lutetium Oxide LU2O3. In 
this case one can cool the nuclear spin subsystem down to very low temperature |l2j and this gives an additional 
enhancement of the effect. In the present work we calculate the coefficient of proportionality between voltage 
and nuclear polarization due to NSM. 

Structure of the paper is the following. In Section II we present simple estimates of the effects induced by 
the electron EDM. Section III contains derivation of the effective interaction induced by the electron EDM. The 
Gd 3+ octupole moment is also discussed here. In Section IV we calculate wave functions of Oxygen electrons 
penetrating inside the Gd 3+ core. Section V presents derivation of the effective Hamiltonian which relates the 
electron EDM with deformation of the lattice. In Sections VI and VII we derive final results for the voltage 
across the sample and the energy shift in the external electric field induced by the electron EDM. Section VIII 
presents similar results for effects induced by the Gd Nuclear Schiff Moment. Finally Sections IX and X contain 
estimates of accuracy of the calculations and our conclusions. 



II. SIMPLE ESTIMATES FOR EFFECTS RELATED TO THE ELECTRON EDM 

In this section we follow simple and important estimates performed in Ref. Q. GdIG is an ionic crystal 
consisting of Gd 3+ , Fe 3+ and 2 ~ ions, see Ref. Uncompensated electronic spins are localized at Gd 3+ 

which has a 4/ 7 electronic configuration and Fe 3+ which has a 4c? 5 electronic configuration. Both these ions 
have half-filled electronic shells, hence the orbital angular momenta of these shells are zero. This explains a 
very weak spin anisotropy in this compound. If the electron EDM is not zero, then both Gd 3+ and Fe 3+ have 
induced EDMs proportional to d e , d a — Kd e . It is known that the coefficient K scales as Z 3 , where Z is the 
nuclear charge, see, e.g. ||. Therefore we neglect the EDM of Fe 3+ and consider only Gd 3+ . GGG differs from 
GdIG in the replacement of magnetic Fe 3+ by nonmagnetic Ga 3+ ions. The most important Gd 3+ ion is exactly 
the same. The enhancement coefficient for Gd 3+ has been calculated previously 

d a = Kd e , If = -4.9 ±1.6 (3) 

If all the uncompensated spins in the sample are 100% polarized then the induced electric field in the sample 
is E = Aimcdda: where ncd = 1.235 x 10 22 cm~ 3 is the number density of Gd. This gives the following induced 
voltage across an L = 10cm sample 

V = EL = Aitn Gd Kd e L = 2 x 10~ 9 V. (4) 

This numerical estimate corresponds to the current limitation on d e , see ([!]) and would be absolutely correct for 
a crystal consisting of neutral atoms or molecules. However the charge of Gd 3+ is not zero. Hence the average 
electric field acting on the ion must be equal to zero, and hence the EDM of the ion cannot really produce the 
field (|J). This means that the crystal lattice relaxes in such a way that the effect (||) is canceled out. This is 
the Schiff theorem Jll| . The theorem is violated due to the finite size of the ion. It means that the estimate (Q) 
should be reduced by some factor depending on the ratio of the ion size to the lattice spacing. Calculation of 
this effect is performed in the following sections. 

On can also perform an estimate of the resulting magnetization due to an applied electric field Eq to the 
sample. The field inside the solid is E = Eq/c where e is the dielectric constant (for GdIG e ~ 15, and for GGG 
e w 30 see Ref. |l|]) The T,P-odd energy shift per Gd ion is given by, 

Se = d a E = Kd e E = 0.9 x 10~ 22 eV. (5) 

The numerical estimate corresponds to the current limitation on d e , see (|l|) and E = lOkV/cm. Using Boltz- 
mann distribution one can deduce the resulting magnetization of the medium. Similar to (Q) the estimate (^) 
contradicts the Schiff theorem. An accurate calculation of 5e is performed below. 



III. INTERACTION OF THE GD 3 + ION WITH ENVIRONMENT, 
VANISHING OF THE OCTUPOLE MOMENT OF THE ION 

In this section we follow the method of calculation suggested in Ref. JTtJ (see also ||). Let 

0(r)= J d 3 R-^L-v-V ext (6) 

be the electric potential produced by the environment of a particular Gd 3+ ion, q(R) is the charge density of 
surrounding ions and E ext is the electric field produced by distant ions and external sources. 
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The energy correction for a system which has a non-zero EDM is 



Ae = ej d 3 r6p(r)<fi(r). (7) 

where Sp(r) is the correction to the charge density of Gd 3+ due to the T,P-odd interaction Htp- 
In equilibrium the average force acting on Gd 3+ vanishes, 

(F Gd3+ ) =-ej d 3 r p V r cj)(r) = 0. (8) 

Here po(r) is the spherically symmetric charge density of Gd 3+ normalized by the condition J ' d 3 r po(r) = 3. 
The total charge density of Gd 3+ is given by p(r) — po(r) + Sp(r). 

Using eqs. (0) and (g) and assuming that E ea;t is uniform in the vicinity of Gd 3+ one can rewrite the energy 
correction as 

Ae = e J d 3 rSp(r)(j)(r) - 1 J d 3 rp (r)(d a • V r )0(r) 

The charge density correction 5p(r) can be expanded in a series of spherical harmonics, 

Sp(r)=Y,fhn(r)Y lm (O r ,<f>r), (10) 

Im 

where 

f m , \- (0\6(R-r)Y lm \n)(n\H TP \0) 

Here we denote the argument of /j m by R to distinguish it from the internal variable r in the matrix elements. 
Because of the negative parity of Htp only odd I contribute to (|Io|). The Y\ m terms correspond to the EDM of 
Gd 3+ , and Y 3m to its octupole moment. First we demonstrate that the octupole and all higher moments vanish 
for Gd 3+ if one neglects the spin-orbit interaction. The kinematic structure of the T,P-odd Hamiltonian is, see, 
e.g- § 

H TP = o(r)s ■ n r , (12) 



where s and r are the spin and the position of the electron, and n r = r/r . An intermediate state |n) in (|ll[) 
belongs to some electronic configuration. We denote the configuration by n. The states within the configuration 
we numerate by an index a, so |n) = \n a ). We also denote by the average energy of the configuration. In 
eq. (|ll]) we first perform summation over a and then over n. With account of this gives 

f fjf\ \^ (0\S(R-r)\n)(n\a(r)\0) 

Ji m {R) oc > [OlYimSinilO) (13) 

_ EjQ — &n 

n 

Here we have used the closure relation ^ Q \n a )(n a \ = 1, where 1 acts in the spin-angular space. This relation 
is valid because of the separation of variables in the problem. Without the spin-orbit interaction, the ground 
state of Gd 3+ reads 

|0) = |L = 0)|5 = 7/2), (14) 

where \L = 0) is the S-wave orbital state (half filled shell). Hence we get from ( fl3| ) 

fim « (0|F JmSi ni|0) = (L = 0\Y lrn n t \L = 0)(S = 7/2\ Sl \S = 7/2). (15) 

It is clear that this equation gives zero for any I higher than 1 because one can combine Yi m with n into total 
angular momentum zero only at I = 1. The spin-orbit interaction 1 ■ s admixes states with nonzero L to the 
ground state of Gd 3+ . The interaction contains the first power of orbital angular momentum 1, so each order 
of the perturbation theory in the interaction changes the selection rule by Al = ±1. Hence, to get a nonzero 
octupole moment of Gd 3+ one has to go to the second order in the spin-orbit interaction. This is a very small 
effect and we do not consider it. 
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Now we return to the energy correction (||). Expanding |r — R| 1 in a series of spherical harmonics 

1 °° 1 4?r r l 

T^J =E E 2|Tl^T5l«(^,0r)y/m(^,^O (16) 



r-. 



Z=0 m=-l ' > 

where r> — max(r, R) and r< = min(r, i?), 



and keeping in mind that according to the previous discussion 5p(r) contains only the first harmonic, we 
transform the energy correction (|^) to the following form 



Ae = e 2 



J d 3 Rq(R)cos6 R i J d 3 r Sp(r)r cos r + R J d 3 



Sp(r) cos 9 r 



(17) 



d 3 rp (r) / d 3 r' 6p{r')r' cos 0' r 



3i? 2 Jo 

This formula is similar to the standard one for the Schiff moment contribution to the T,P-odd energy correction, 
see Refs. [^7|J^1. For a smooth function q(R) the eq. ( |l7T) can be easily transformed to the standard form. However 
the point is that in the present case the density q(R) is not a smooth function of R. Similar to the standard 
case, eq. (0) gives zero Ae if the external charge q(R) is localized outside of the Gd 3+ core. So we have to 
calculate the probability for external electrons to penetrate inside the Gd 3+ core. This calculation is performed 
in the following section. Considering the integral over R in eq. (17), we deduce that only an asymmetric part 
in q(R) proportional to cos#a. can contribute to the energy correction. There is no such term for the lattice 
without deformation. Hence the nonzero Ae is also related to the deformation of the lattice. 



IV. PENETRATION OF OXYGEN ELECTRONS INSIDE THE GD 3+ CORE, DEFORMATION OF 

THE LATTICE. 

In the Garnet structure any Gd 3+ ion is surrounded by eight 2 ~ ions organized in a dodecahedron structure 
(slightly distorted cube) p3[ . For the calculation we approximate this structure by a cube with Gd 3+ in 
the middle. The Gd-0 distance used is 4.53as, where as is the Bohr radius. In the case of GdIG there 
are also uncompensated spins localized at Fe 3+ ions. However the T,P-odd effect scales as Z 3 , and therefore 
the contribution of Fe is negligible compared to that of Gd. An additional reason for suppression of the Fe 
contribution is that the crystal radius of Fe 3+ is much smaller than that of Gd 3+ . The quantity of interest 
is the charge density of external electrons q(r) inside the Gd 3+ ion. The electrons are coming from the eight 
nearest 2 ~ ions. 

The O 2- ion has a closed shell configuration, with 3 double occupied 2p-orbitals. One can show that 2p v 
orbitals do not contribute to the effects we are interested in. So we only consider 2p CT -orbitals pointing towards 
Gd 3+ as schematically shown in the 2-dimensional picture in Fig. [|. Hence there are 16 electrons that can 
potentially penetrate inside Gd 3+ . 




FIG. 1. A schematic 2-dimensional picture for penetration of 2p a -orbitals of O 2 inside Gd 3 
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First of all let us rewrite the eight 2j3 CT -orbitals in the basis with definite symmetry with respect to the cubic 
structure. The new basis reads 
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= (1/V8)( 
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Here, for example |1) denotes the 2p CT -orbital of the first oxygen. The ions are numbered as in Fig. g. All 
of the states (|18| ) are doubly occupied. As long as the electron is moving near the Oxygen ion the electron 
wave function can be described as a 2p (T -orbital. However when the electron approaches the Gd 3+ ion its wave 
function is strongly deformed by the potential of Gd 3+ . To describe this effect we draw an imaginary sphere 
of radius T around the Gd 3+ ion as shown in Fig. jjl W e assume that outside the sphere the Oxygen potential 
dominates and the wave is described as one of the lfrs|) orbitals. However inside the sphere the Gd 3+ potential 
dominates and hence the same orbital is described as s-, p-, d-, or f-orbital of Gd. We are interested only in s- 
and p- orbitals that correspond to \S), \P X ), \Py), and \P Z ) states from (filf). At r < Y we will use the 6s and 
6p orbitals of Gd 1+ as the basis. We can also choose 7s and 7p or 8s and 8p, the final result is independent of 
the choice because at r < Y all these wave functions are proportional to each other. 




FIG. 2. The matching sphere around Gd with the surrounding eight Oxygen neighbors. 



Thus at r < T, we get 

|S)^/?s|6 S }, \P Z ) -» P \6p z ), \P X ) - Pp\6 Pz ), \P y )^(3p\6 Pz ). (19) 

To find the coefficients fis and ftp we match (|l^) with ([l8]) at the surface of the sphere with radius r — Y using 
the following projection 

(u(T,d),e)\v(T,d),e)) r = / d<t> / d6 sinO u*(T,(j>, 6) v{T,<f>, 6). (20) 

Hence 

_ (6s|S) r _ (QPz\P)r m , 

/S <6s|6 S ) r PP (6p 2 |6p z )/ 1 1 

This is a simple and reliable method to calculate the coefficients (3s and j3p. It has been previously used for 
calculation of the nuclear quadrupole resonance frequencies in La2Cu04, see Ref. Jl8| ]. In that case the accuracy 
of the method was verified by experimental data and proved to be ~10%. For numerical calculations with 
(p(i|) , ( 2l|) we have used 2p- wave- functions of 2 ~ calculated previously in The 6s and 6p wave functions 
of Gd" have been obtained by the relativistic Hartree-Fock method for configuration Is 2 ...Af 7 6s6p. In this 
calculation we average over polarizations of the unclosed shells. It means that 4/ 7 is taken as 4/ 14 with 
50% population of each orbital, similarly 6s is taken as 6s 2 with 50% population and 6p is taken as 6p 6 with 
1/6=16.6% population. There is also the fine structure effect: the difference between 6pi/2 and 6^3/2 states. 
However this effect is small and has been neglected in the calculation of (3p. The values of (3s and (3p calculated 
with (po|),(|2l]) for different choices of the matching radius Y are presented in Table 1. The coefficients are not 
sensitive to the value of Y and this confirms validity of the method. 
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TABLE I. The coefficients /3g, /3' s , /3p, and j3' P for expansion of Oxygen 2p-orbitals in terms of Gd 1+ 6s- and 
6p- wave- functions, see eqs. (p7j),(|l9|). The coefficients are presented for different values of the radius of the matching 
sphere F. P s and P p represent the probability for the corresponding Oxygen electron to penetrate inside the matching 
sphere. (N2 P ) is the average number of Oxygen 2p-electrons inside the matching sphere. 

The calculation assumes that the 2p-wave-functions of O 2- are "rigid" , i.e. at r > T they are not influenced 
by Gd. This can be true only if the total probability for an Oxygen electron to penetrate inside the matching 
sphere, r < T, is small. To check this we calculate the probabilities 



Ps = Ps / d*r\^ s {r)\\ 



P p =0 2 p [ d 3 r\ij 6p (r)\ 2 
Jo 



(22) 
(23) 



Numerical values of P s and P p are presented in Table 1, they are really very small. The average number of 
Oxygen electrons penetrating inside the matching sphere, 



(N 2p ) =2(P S +3P P ), 

listed in Table 1 is also rather small. Thus, the external charge density with respect to Gd 3 " 
sphere is given by 

g(r) = <zo(r) = 2/?^ 6s ( r )| 2 + 2 # (|V> 6p ,(r)| 2 + \^ Pv (r)\ 2 + |^ 6p2 (r)| 2 ) 



(24) 

inside the matching 
(25) 



The coefficient 2 is due to the double occupancy (spin up, down) of each orbital. 

It has been pointed out in the previous section that the charge density (^5|) is insufficient to generate a T,P- 
odd effect because <?o( r ) is symmetric with respect to reflection r — * r. To have an asymmetric part Sq(r) one 
needs to consider a deformation of the lattice. We will assume that Gd is shifted along the z-axis with respect 
to the environment by a small displacement X, for definition of axes see Fig. ^. The displacement is shown 
schematically in the 2-dimensional picture in Fig. [|. In the general case the displacement can have an arbitrary 
direction. However we do not expect that the final answer will be very sensitive to the direction. Therefore we 
consider the simplest geometry: the displacement along the axis of the cube. 









flu/ / 

If 


^4 





FIG. 3. A schematic 2-dimensional picture for penetration of 2p a -orbitals of O 2 inside the shifted Gd 3+ 
It is clear that outside the matching sphere nothing is changed: eqs. (|8[) describe the wave functions. 
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However inside the sphere the wave functions are changed compared to (tL9j 



\S)^p s \6s) 



as 



\P Z ) -> I3p\6p z 



I3' P — |6s) 
a B 



(26) 



We present only the orbitals that give a nonzero contribution to the T,P-odd effect. To have dimcnsionlcss 
coefficients (3' s and j3' P we divide the displacement by the Bohr radius. To calculate the coefficients f5' s and (3' P 
we use the same matching procedure (E^ 



d (6^|5) r 



dX (6 Pz \6 Pz 



0'p 



d (6s\P z 



dX (6s|6s) 



(27) 



The values of j3' s and /3p for different values of the radius of the matching sphere T are presented in Table 1. 

As one can see from Table 1 the values for (3g, (3p, j3' s and j3' P are almost independent of the choice T in 
the interval 2.2as — 2.8as, this is the dual description region. Existence of this region confirms validity of 
our approach. For all further calculations we will use T = 2.han- Using eqs. (^6|) we immediately obtain the 
asymmetric term for the charge density, 



8q(r) = 4(PsP's 



PpP'p) — -0 6s (r)-06 Pz (r) 
ap 



(28) 



In the above calculation we have neglected relativistic effects. This is justified as far as the calculation of the 
coefficients (3 is concerned. However for the behavior of Sq(r) near the Gd nucleus the relativistic effects are 
very important. Therefore we have to modify ( f?8|) to account for this effect. Fortunately this modification is 
obvious and immediately follows from decomposition of 6p z in terms of 6pi/2 and 6p3/2- This gives 



5q(r) -> 4(Ps(3' s + PpP'p) 




P3/2 



Pl/2 



(29) 



where 



X Pl/2 



-^P 1/2 1/2 



X P3/2 



Rj>3/2 ^P3/2 
*^P3/2 ^P3/2 



R s il s 

iR s Q s 



(30) 



are the Dirac wave functions of 6pi/2, &P3/21 and 6S1/2 states. f2 and fi = — (cr • n)Q are the usual two-component 



spherical spinors |19|; R and R are upper and lower radial wave functions that have been calculated using the 
Relativistic Hartree-Fock method. 



V. RELATION BETWEEN THE LATTICE DEFORMATION AND THE T,P-ODD ENERGY 
CORRECTION INDUCED BY THE ELECTRON EDM 



The T,P-odd interaction of the electron EDM with the electric field of the Gd nucleus E is Vd — — d e 7o£ ■ E, 
and account of the Schiff theorem for electronic degrees of freedom reduces it to the form, see e.g. Ref. B 



Vd -> VJ = -d e ( 7o - 1)33 • E, 



(31) 



where 70 and £ = 70757 are Dirac 7-matrices. According to [Q the EDM of Gd 3+ ion is saturated by 4f-5d 
mixing 



d a = Kd e = e / d 3 r 5p{v)r cos r w 2e 

J m=-: 

Diagrammatically it is shown in Fig. 



(4/ m I r cos r 1 5d m ) (5d m \ Vf\4f m ) 
E&f — E 5( i 



(32) 



4f 



5d . 4f 



4f 5d 



4f 



FIG. 4. Leading contribution to Gd 3+ EDM. The dashed line denotes the dipole moment ez = ercosO, and the cross 
denotes the T,P-odd interaction VJ , eq fclA)- 
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The same 4f-5d virtual transitions saturate the T,P-odd correction to energy (17), and taking into account eq. 
( ^9| ) the correction (0) can be rewritten as 



a 4e 2 , , X * (r 3 m|co S g|y4 m )(5d m |^|4/ m/ , 
A e = ^(PSPs + l jp PP)~ 2 2^ ET~ — (Ji + h + h) 



h + h) ■ 



(33) 



where, 




h = 


S? G 


h = 


/ dRR 3 
Jo 


h = 


i r r 


( r fd) 


= / dr' r 
Jo 



u Pl/2 



RsRp 1/2 ) + ^(RsRp 3 / 2 



L Pl/2 



L Pl/2 



L Pl/2 



L Pl/2 



■i- 



R S R 



/2 . 



L Pa/2 



L P3/2 



^P3/2 RsRp 3 /2 ) 



drr 3 R 4f {r)R 5d (r), (34) 



dr R4f(r)R 5 d(r), 

R r r 

d 3 rp (r) / dr' r' 3 R Af {r')R 5d {r'), 
o Jo 



In the integrals we have replaced the upper limit of integration by Y = 2.5as, because eq. (|29|) is valid only 
at r < r. The final result is not sensitive to the upper limit because these integrals are convergent at smaller 
distances. The numerical values obtained with Hartree-Fock wave functions are 



h = -4.1 x 1CT 3 «„' 



h = -5.3 x 1CT 3 a B \ h = -7.5 x 1CT 3 a B \ (r fd ) = 0.55a B . (35) 



Substituting these values, as well as the value of the coefficient K from (||) and the coefficients (3 from Table 1. 
in eq. ( |33| ) we find the T,P-odd energy correction per Gd ion as a function of the ion displacement with respect 
to surrounding Oxygen ions. 



Ae(X) 



e X 

~ J\ Ty Cfcg 

cig a B 



S a% 



with A = 0.11 . 



(36) 



The vector form is valid for an arbitrary direction of the displacement X and an arbitrary direction of the 
Gd 3+ polarization S/S. The integrals ( pq ) give A = 0.09, however they only account for the contribution of 
2p-orbitals of Oxygen. We know jl8j that there is also a contribution of 2s-orbitals of Oxygen that is 20%-25% 
of the 2p-contribution. This is why we take A = 0.11. 

The integral 1% corresponds to the compensating term in (jl7|) and the contributions corresponding to the 
integrals l\ and I2 can be represented diagrammatically as shown in Fig.[| 



6s 
4f 



5d 



6p 
4f 



+ permutations 



FIG. 5. Direct contribution to the T,P-odd energy shift. The dotted line denotes the Coulomb electron- electron inter- 
action, and the cross denotes the T,P-odd interaction VJ , eq 



This is what we have taken into account and these are direct diagrams only. In principle there are also exchange 
diagrams shown schematically in Fig. ^| that we have neglected. 



6s 
4f 



5d 



4f 



6p 



+ 



4f 
6s 



6p 



+ permutations 



FIG. 6. Exchange contribution to the T,P-odd energy shift. The dotted line denotes the Coulomb electron- electron 
interaction, and the cross denotes the T,P-odd interaction VJ , eg d?H). 



The exchange diagrams have higher angular momenta of the Coulomb quantum than the direct ones, the 
diagrams in Fig|] have k = 2,3,4, while the direct diagram in Fig jE] has k = 1. Therefore it is unlikely 
that the exchange contributions can be comparable with direct ones. Nevertheless, at some later stage, these 
contributions have to be calculated as well. 



VI. T,P-ODD VOLTAGE ACROSS A MAGNETICALLY POLARIZED SAMPLE OF GDIG 

In this section we consider a GdIG sample magnetically polarized along some axis. Then according to eq. 
( |36| ) each Gd ion can gain energy from a small distortion of the lattice. However the lattice has stiffness and 
therefore the total variation of energy per Gd ion as a function of the displacement X is of the form 

Ae T = \k* X 2 - A^—d e , (37) 
l a B a B 

where K e i is the effective elastic constant per Gd ion. Minimizing (|37]) with respect to the lattice deformation 
X we find the deformation induced by the electron EDM. 

X = -^rde. (38) 

To calculate the deformation we need to know the value of K e i. The elastic constant is related to optic phonons, 
and can be found from the spectra of the phonons. However we will use a simpler method based on the known 
static dielectric constant of GdIG, e w 15 jl(|. Consider some external electric field applied to the GdIG sample. 
According to standard relations 

D = eE = E + 4ttP, (39) 

where E is the average electric field in the sample and P the polarization per unit volume. Assuming simple 
cubic structure the local field acting on each ion is [^0| 

E 1 = E+^P = 47rP + (40) 

This gives the displacement of the Gd 3+ ion in the field 

3eEi Uire „ / 1 1\ 

Similarly one can find the displacement of the Fe 3+ ion X2 in terms of the corresponding elastic constant K' e[ . 
On the other hand, the polarization per unit volume by definition is P — Zexinod + 5ea;2«Gd- Altogether this 
gives the following result for K e i. 

We recall that nod = 1.235 x 10 22 cm~ 3 is the number density of Gd in GdIG. For further estimates we will 
assume that K e i = K' el . There is no reason for these constants to differ substantially, and moreover eq. ([42| ) is 
not that sensitive to the ratio K e i/K' el . 

Coming back to the situation without any external electric field and using eqs. (|3^) and (f42| ) we find the 
electric polarization P = ieXncd and hence the electric field inside the sample 

A d 

E = -AttP = r • -I- = -0.54 x 10- 10 F/cm. (43) 

H^t + ^J aB 

This numerical value corresponds to the current upper limit on the electron EDM (|l|). For a 10 cm sample it 
gives a voltage AV = 0.54 x lO -9 ^. This is smaller than the naive estimate (||) and this is suppression due to 
the Schiff theorem. Fortunately this suppression is not strong. 
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VII. T,P-ODD MAGNETIZATION OF GGG IN THE EXTERNAL ELECTRIC FIELD 



Following the suggestion of Lamoreaux |l| we consider now a different kind of experiment: an external electric 
field is applied to the sample and this leads to the macroscopic magnetization. Equation ( [4l| ) gives the Gd shift 
due to the electric field. 

X = Jf ~ X) E (44) 
327rencd 

Substituting this into eq. (|3|) we find the energy shift of the Gd 3+ ion 

Ae = - w(g_g) _ 2 . 8 x 10 -- e y. (45) 

32n(a%n G d) S 

The numerical estimate corresponds to E — lOkV/cm, the current limitation on d e , ([!]), and the dielectric 
constant e s» 30 for GGG Jl6| . ft is interesting that in spite of the Schiff theorem the result (^) is larger than 
the naive estimate ( a). 

The energy shift ( |45| ) leads to the macroscopic magnetization of the sample. The magnetization depends 
also on temperature and internal magnetic interactions in the compound. We do not discuss these points here. 
According to estimates |jj the magnetization due to the energy shift ~ lO~ 22 eV can be measured and moreover 
the prospects for improvement of sensitivity are very good. 



VIII. T,P-ODD VOLTAGE AND MAGNETIZATION DUE TO THE NUCLEAR SCHIFF MOMENT 

The effects considered above are due to the electron EDM, so they are sensitive to T,P-violation in the lepton 
sector. Similar effects are generated by the Nuclear Schiff Moment (NSM), so they are also sensitive to T,P- 
violation in the hadron sector. We begin from consideration of the same compounds GdlG and GGG, assuming 
only one Gd isotope. This must be an odd isotope, say 155 Gd. The NSM is a T,P-odd vector moment that 
is a property of the nucleus, see Refs. p7|,p| 

= S N j, (46) 

where I is nuclear spin. In previous sections we denoted by S the total electron spin, note that it has nothing 
in common with the NSM Sjv- Due to the Schiff moment there is the T,P-odd interaction between electron and 
nucleus. This is the contact interaction, therefore the matrix elements of the interaction are nonzero only for s- 
and p-electrons. These matrix elements have been calculated in Ref. |r7| , see also Q 

< s l^> = 4 f efl wa S*r(n,|nin p ). (47) 

Here Q s and Q p are spherical spinors corresponding to the s- and p-states; n is the unit vector; v s and v v are 
principal effective quantum numbers; and R is the relativistic factor that is different for pi/% and p 3/ / 2 electrons 

a 271/2-2 

_ 47i/2^ . , 

^1/2 ~ tftk:, — r-rvi2> ( 48 ) 



#3/2 



[r(2 7l / 2 + i)] 2 ' 

48 7l /2^ 1/2+73/2 ' 3 

r(2 7l/2 + i)r(2 73/ 2 + 1) ' 



Here 7l / 2 = Vl — Z 2 a 2 , 73/2 = v4 — Z 2 a 2 , a is the fine structure constant; T(x) is the usual 7 -function; and 
xq = (2Zro/a,B), where r$ is the nuclear radius. For Gd, R1/2 = 3.56 and -R3/2 = 2.99. Assuming that the 
nucleus is completely polarized along the z-axis, and taking also \p) — ► \p z ) we find from ( [47|) 

According to (^6|) lattice deformation X leads to the mixing of 6s- and 6p-orbitals. Hence the matrix element 
of Hs over the electronic states of the deformed lattice is equal to 

(H S ) = * (fo^+foflO— u NZ \ /2 (\Ri/2 + 1^/2] ■ (50) 
7rV3 a B a%{v s v p yi 2 \3 1 3 ' 
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The effective principal numbers for 6s- and 6p-orbitals are: vq s = 1.24, ^6 Pl/2 = 1-46, and vep 3 , 2 = 1-49, see 
comment pi] ]. Substituting these values, together with coefficients (3 from Table 1 and all the other parameters 
we find the T,P-odd energy correction due to the NSM 

Ae(X) = -B-r- — S N -> -B-^(X-I) with B = 1,2 x 10 3 . (51) 
a% a B I a B 

The vector form is valid for an arbitrary direction of the displacement X and an arbitrary direction of the Gd 
nuclear spin I. Equation ( |5l| ) for the energy shift due to NSM is similar to eq. for the energy shift due to 
the electron EDM. We stress once more that Sn in denotes the NSM while S in (|3^) denotes the electron 
spin. In (51) we have also taken into account 20% contribution from the Oxygen 2s-electrons, see comment 
after eq. (3£). 

Next we follow the same path as we did for the effects due to the electron EDM. First consider a sample with 
fully polarized nuclear spins, no external electric field. Then, repeating the calculations performed in Section 
VI we find the following electric field induced inside the solid 

E = ^ • ^- = -0.9 x \Q- 12 V/cm. (52) 

~ 3 

The numerical value corresponds to the current upper limit on the NSM of 199 Hg (||). 

If the external electric field is applied to the solid then repeating the calculations from Section VII we find 
the energy shift per Gd depending on the orientation of Gd nucleus 

Ae = - ^~ l \ • 4 • - -0.47 x W^eV. (53) 

327r(a|n G d) a% I 

The numerical estimate corresponds to E = lOkV/cm, the current limitation on NSM of 199 Hg (|^), and the 
dielectric constant e ss 30 for GGG llql . 



IX. ACCURACY OF THE CALCULATIONS 



The effects related to the Nuclear Schiff moment are sensitive only to s- and p-orbitals of Gd. There are 
three sources of possible corrections to the results presented in eqs. (|5l|), (|52"|), and j5^): 1) interatomic many- 
body correlations, 2)inaccuracy of the used matching procedure for calculation of the expansion coefficients 
P, 3)deviation of the effective elastic constant K e i for Gd from that for Fe. Usually for s- and p-orbitals the 
correction due to the interatomic correlations does not exceed 20%, see e.g. [||. As we already mentioned above 
the matching procedure has been previously used for calculation of the nuclear quadrupole resonance frequencies 
in La2Cu04, see Ref. iQ. In that case the accuracy of the procedure was verified by experimental data and 
proved to be about 10%. We believe that this estimate is also valid for the present case. Finally, even in the 
unlikely case of 40% difference in elastic constants K e i and K', the variation in (|5^) and (|5^) is 20%. All in all 
this gives a 20%-30% estimate for a possible inaccuracy in (f5l[), (|52]), and (|53|). 

Situation with effects induced by the electron EDM is more complex. The mechanisms 1-3 for possible 
corrections contribute here as well. The interatomic many-body correlations have been considered in Ref. ]l4| . 
They contribute at the level 20%, but the most important ones have already been accounted in the value of 
the EDM enhancement coefficient K we used in the present work. In addition there are uncertainties due 
to 4)experimentally unknown E±f — E$d energy splitting in Gd 3+ , and 5)unaccounted exchange diagrams in 
Fig. 6. The most important is the energy splitting E^f — E^ that determines uncertainty in the electron EDM 
enhancement coefficient K, see eq.(||). We estimate the accuracy of the results (p6|), ([i"3|), and ( |45| ) as 30-40%. 
Experimental and/or theoretical determination of the £4/ — E^ energy splitting in Gd 3+ would be the most 
important to improve the accuracy of the calculation of the effects related to the electron EDM. 



X. CONCLUSIONS 



In the present work we have considered the T,P-odd effects in solids. The calculations are performed for 
Gadolinium Gallium Garnet, GdaGasO]^ (GGG), and Gadolinium Iron Garnet GdaFesO]^ (GdIG). We consider 
the effects due to the electron electric dipole moment (EDM) and due to the Nuclear Schiff Moment (NSM). 
Both GdIG and GGG have uncompensated electron spins on Gd 3+ ions. There are two possibilities to probe the 
electron EDM. The first one is to polarize magnetically the electron spins and to measure the induced voltage 
across the sample. According to our calculations at the current limitation on the electron EDM, (R), the induced 
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voltage across a 10cm sample is 0.54 x 10 9 V, see eq. (|3|). Another possibility is to apply an electric field to 
the unpolarized sample. This leads to the spin-dependent energy shift of each Gd ion Ae = 2.8 x 10" 22 el/ at 
E = lOkV/cm, see eq. ([l5"|). This can be measured via macroscopic magnetization of the sample. 

Gd nuclei can be polarized by the electron polarization via the hyperfinc interaction or they can be polarized 
independently. For 100% polarization at the current limitation on the NSM (^) the induced voltage across a 
10cm sample is 0.9 x 10 _11 V^, see eq. (|52]). Application of the external electric field to the unpolarized sample 
leads to the energy shift of each Gd Ae = 0.47 x 10~ 23 eV at E = WkV/cm, see eq. (|||). The shift is dependent 
on the nuclear spin. This can be measured in NMR experiments or via macroscopic electronic magnetization 
of the sample due to the hyperfine interaction. 

Another possibility for searches of NSM is to use compounds without unpaired electrons, such as Lutetium 
Gallium Garnet, Lu3Ga 5 0i2 or Lutetium Oxide LU2O3. In this case both the NMR and the macroscopic 
magnetization techniques can be used. A possibility to cool the nuclear spin subsystem down to very low 
temperatures 12 can give an additional advantage. 

In our estimates of the effects induced by NSM we have used the current limitation (j|) for NSM of 199 Hg 
as the reference point. However the NSM is dependent on nucleus. 199 Hg is a simple spherical nucleus while 
155 Gd is a deformed nucleus with close (105keV) levels of opposite parity and low energy (~ lMeV) collective 
3 _ excitation. In this situation an additional enhancement of NSM by an order of magnitude is possible. This 
problem requires separate consideration. 

There is one more mechanism to probe T,P-violation inside nucleus. It is related to the T,P-odd nuclear 
electric octupole moment and it also requires separate consideration. This mechanism was mentioned in Ref 
p7[ , but has never been considered. For 155 Gd it can give a large contribution because of the low energy 
collective 3~ excitation. 



ACKNOWLEDGMENTS 



We are grateful to S. K. Lamoreaux, L. R. Hunter, and D. Budker who attracted our attention to the problem. 
We are also grateful to them as well as to W. R. Johnson and M. G. Kozlov for stimulating discussions and 
interest in the work. 



[1] S. K. Lamoreaux, |nucl-ex/0109014 

[2] L. R. Hunter, T alk at the workshop Tests of Fundamental Symmetri es in Atoms and Molecules, Harvard, 2001. 

Available online http://itamp.harvard.edu/fundamentalworkshop.html 
[3] J. H. Christensen, J. W. Cronin, V. L. Fitch, and R. Turlay, Phys. Rev. Lett. 13, 138 (1964). 
[4] A. D. Sakharov, Sov. Phys. JETP Lett. 5, 24 (1967). 

[5] J. Schwinger, Phys. Rev. 91, 720, 723 (1953); W. Pauli in Niels Bohr and the Development of Physics (Pergamon 
Press, London, 1955). 

[6] I. B. Khriplovich and S. K. Lamoreaux, CP Violation Without Strangeness (Springer, Berlin, 1997). 

[7] B. C. Regan, E. D. Commins, C. J. Schmidt, and D. DeMille. Phys. Rev. Lett. 88, 071805 (2002). 

[8] D. DeMille. Tal k at the workshop Tests of Fundamental Symmetrie s in Atoms and Molecules, Harvard, 2001. 
Available online http://itamp.harvard.edu/fundamentalworkshop.html 

[9] F. L. Shapiro, Usp. Fiz. Nauk 95, 145 (1968) [Sov. Phys. Usp. 11, 345 (1968)]. 
[10] B. V. Vasil'ev and E. V. Kolycheva, ZhETF 74, 466 (1978) [Sov. Phys. JETP 47, 243 (1978)]. 
[11] M. V. Romalis, W. C. Griffith, J. P. Jacobs, and E. N. Fortson, Phys. Rev. Lett. 86, 2505 (2001). 
[12] A. S. Oja and O. V. Lounasmaa, Rev. Mod. Phys. 69, 1 (1997). 
[13] Physics of Magnetic Garnets, ed A. Paoletti, Nort h-Holland, Amste rdam, 1978. 
[14] S. Y. Buhmann, V. A. Dzuba, and O. P. Sushkov, |physics/0204076 . 
[15] L. I. Schiff, Phys. Rev. 132, 2194 (1963). 

[16] The value e ~ 15 for GdIG is based on data for Y-Gd Iron Garnets retrieved March 1, 2002, from 



components.com/temex/catalog/ciron.pg) 

The value e w 30 for GGG is based 



http: / /www.temex- 



on data retrieved March 1, 2002, from http://www.mt 



[17] 

[18] 
[19] 
[20] 



berlin.eom/charts/chart_07.htm#IND7 

O. P. Sushkov, V. V. Flambaum, and I. B. Khriplovich, Zh. Eksp. Teor. Fiz. 87, 1521 (1984) [Sov. Phys. JETP 60, 
873 (1984)] 

V. V. Flambaum and O. P. Sushkov, Physica C 168, 565, (1990). 

V. B. Berestetskii, E. M. Lifshitz, and L. P. Pitaevskii, Relativistic quantum theory (Pergamon Press, Oxford, 1982). 
C. Kittel, Introduction to Solid State Physics (John Wiley & Sons, New York, 1971). 



12 



[21] We found the effective principle numbers by direct comparison of the Relativistic Hartree-Fock wave functions at 
small distances with the corresponding semiclassical formulas for a neutral atom, see, e.g. Ref. |J. This is why we 
do not account for the effective charge of the ion Zi = 2 in eqs. (^),(^), and (^). One can account for Zi, but 
then the effective principle numbers have to be rescaled, v — »■ fx: l/v 3 — Zf /fj?. The matrix element is invariant 
under this procedure. 



13 



